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Abstract 

We develop the hnear equations that describe graviton perturbations of ArfS's-Schwarzschild 
generated by a string traihng behind an external quark moving with constant velocity. Solv- 
ing these equations allows us to evaluate the stress tensor in the boundary gauge theory. 
Components of the stress tensor exhibit directional structures in Fourier space at both large 
and small momentum. We comment on the possible relevance of our results to relativistic 
heavy ion collisions. 
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1 Introduction 

If an external quark (meaning an infinitely massive, fundamentally charged, point-like par- 
ticle) is passed through a thermal plasma of SU {N) A/" = 4 super- Yang- Mills with both 
and Oym-^ large, it experiences a drag force [TJ |2] 

dt 2 v/TT^ ' ^ ' 

where v is the speed of the quark and T is the temperature of the plasma. The result 
was derived using the duality between string theory in anti-de Sitter space and conformal 
field theory on the boundary of anti-de Sitter space (AdS/CFT) IH E]- The diffusion 
constant D = 2/ {-nT gy implied by (H)) was derived independently in 0, also using 
AdS/CFT. A number of other papers 0131111011111111113111111311211131111111201 use 
AdS / CFT in related ways to describe energy dissipation from a moving quark. This work is 
motivated in part by the phenomenon of jet-quenching observed in RHIC experiments: for 
recent experimental accounts see for example [211 1221 12S] • 

It was observed in |2] that (T^n) in the boundary gauge theory is calculable via 
AdS / CFT but requires a technically non-trivial analysis of graviton perturbations in AdS^- 
Schwarzschild. The aim of this paper is to develop the relevant equations and solve them. 
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both in limits that are analytically tractable and through use of numerics for selected values 
of the velocity. Section |21 comprises a derivation of the equations. The linearized gravi- 
ton equations are stated in full at the end of section |21 Section El includes solutions of the 
equations near the boundary and near the horizon of AfiS's-Schwarzschild, as well as ex- 
pressions for the stress tensor in the near-field limit. A boundary- value problem in classical 
five-dimensional gravity is stated at the end of section 13.21 which determines (Tm„) in the 
boundary theory. Section 0] presents the results of numerical work, and section El is devoted 
to discussion of the possible relevance to jet-quenching in relativistic heavy-ion collisions. 

The reader wishing to skip all the five-dimensional technicalities and see the "answers" 
may skip directly to section jH with the following definitions and conventions in mind: 

• Wave- numbers K = [Ki, K2, K3) are rendered dimensionless by including a factor of 
Zh = l/vrT. 

• Usually we set K3 = and K2 = K± > 0. 

• Often we refer to K = ^Kf + Kl and 9 = tan~\K±/Ki). 

• C^mn) is the i^-th Fourier coefficient of the co-moving part of {Tmn) in the 3-1-1- 
dimensional boundary gauge theory: see (jl^ . 

• is a dimensionless quantity proportional to (T^^): see 

• Using symmetries and conservation laws, the non-zero components of Q^^^ can be 
expressed in terms of three complex quantities Q^*, Q^ff', and Q^*: see (jHTj) . 

• is the easiest quantity to interpret, as it is directly proportional to (T^), the K-th 
Fourier component of the energy density. 

• Qx for X = A, D , OT E is Qx* with the Coulombic near-field of the quark subtracted 
away: see ()62|). Note also that the inhomogeneous term Pmn in amounts to far- field 
subtraction in the definition of both and Qx- 

2 The graviton equations of motion 

The relevant part of the action for supergravity plus the string is 




X 



V^{R+12/L^) 
24 



27ra' 



1 



/ 



(2) 
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where L is the radius of AdS^. We have excluded the dilaton from (j2I) because it decou- 
ples from the metric at the level of linear perturbations around the AdS's-Schwarzschild 
background, 

dsjo) = G^l^ldx^dx"" = {-hde + dx" + h = l-y\ (3) 

We have introduced in (jHl) a radial variable y which runs from at the boundary of AdS^- 
Schwarzschild to 1 at the horizon. A more conventional choice of radial variable is z = zhV- 
The position of the string can be described in static gauge as 



Xf^it,y)= (t X\t,y) y 
X\t, y)=vt + ^y) ay) = flog ^ + ^ log ^ + 

The equation of motion is 



1-y 



(4) 



R^"" - -G^'^R - -^G^" = T^"" , (5) 



where 

.2 



2 p I 

27ra' J 



..2 3 
^5 



(6) 



--^y^Vl^5{x^ -vt- ay))6{x^)6{x^)daX''d''X' 



2'Ka' L3 " 

is the stress-energy tensor due to the string, expressed explicitly in static gauge in the second 
line. If we perturb 

G,, = Gf} + h,,, (7) 
then, schematically, the form of the linearized equations following from is 

^Adsh^"" = r^" (8) 

where h^^ = G^^^-^G^^^-^hp^ and A^^^ is a variant of the Lichnerowicz operator. 

The stress tensor r^*^ depends on x^ and t only through the combination x^ — vt. Thus 
we can expand 

T^\t,x\x\x\y) = [ f^ r^-^y)^[J<.i-'-t)+K,.^+K,.^]hu (9) 
J [2tx) 
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(Note that K = znk = k/nT is dimensionless.) If our interest is the co-moving graviton 
response, we can make a similar expansion for h^^. Then, again schematically, one obtains 
from (jHJ a set of coupled ordinary differential equations in y: 



0. 



(10) 



It is these equations which we wish to formulate more precisely and then solve. Note that 
although the Fourier modes and h'^ are complex, they satisfy conditions like r^'^ = 
{r^y that ensure the position space quantities are real. From the asymptotic behavior of 
h'^ near the boundary of Aci^'s-Schwarzschild one may extract the A'-th Fourier mode (T^„) 
of the co-moving contribution to the stress tensor. A detailed discussion of the extraction of 
{T^^) is deferred to section ITTl 

The rotational symmetry around the axis of motion of the quark enables us to choose 
K = {Ki, K±,0) with K± > 0. The remaining symmetry is a Z2 sending — —x^ 
metric perturbation can be parametrized as 



The 



''liu 







H02 


-^^03 


o\ 




Hn 


H12 


-^13 





H20 


H21 


H22 


H23 







H31 


H32 


H33 





\ 











0/ 



The vanishing entries represent a gauge choice which we will refer to as axial gauge. The 
J^-th Fourier mode of the string's stress tensor is 
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(,-iKili(y)/z„ y5 
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(12) 



where ^{y) is as given in 

Because (llUj) is an equation for symmetric rank-two tensors, it has 15 algebraically in- 
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dependent component equations. Ten of these, namely = with < m < n < 3, are 
second order, and the other five, S'^^ = with fi unrestricted, are first order constraints.^ 
There are 10 dependent variables Hmn, so the full system S'^'^ = seems over determined. 
But it isn't: if the constraints are imposed at one value of y and the second order equations 
are then solved, the constraints continue to hold automatically for all y. 

The Z2 symmetry that takes — > —x^ causes the equations (jiup to partially decouple. 
The Z2 "charge" of a component Hmn of the metric perturbation is the parity of the number 
of indices equal to 3. Similar charge assignments can be made to the equations S'^'^ = 0. 
Thus for example Hq^ and S^^ = are odd while and S^^ = are even. The three 
odd second order equations of motion and the one odd constraint equation involve only the 
odd variables, whereas the seven even second order equations and the four even constraint 
equations involve only the even variables. Moreover, only the even equations involve non- 
zero components of the stress tensor (fT^ . So it is consistent to set all the odd variables 
equal to from the outset. In the interests of generality, we will not do this yet, but rather 
consider how the equations may be further decoupled. 

Briefiy, we make the following definitions and find the following differential equations: 

K = JkI + KI e = tan-i (13) 



A 



~Hii + 2 cot eHi2 - cot2 eH22 + csc2 OH^s 
2^^ 



(14) 



h 



(15) 



03 



Bo 



Hrs + tan ^^23 



(16) 



+ f 




-h v'^cos^ 9h 
-1 v'^ cos^ 6 




(17) 



B[ -hB'^ = (18) 

^Wc associate y with /j, = 5. There is no coordinate associated with /i = 4. This convention serves as a 
reminder that y is the fifth dimension. 



5 



c 



sin + cos Qil2?, 
K 



y h J 



C = 



Hoi-cot9Ho2 n -Hu + 2cot29Hi2 + H22 

J^l — 7. J^2 — 



2v 



(19) 
(20) 
(21) 



y 



-h f^cos^^/i 



^M-1 v'^cos^e 



D2. 



D[ - hD'2 



ivKi 



(22) 
(23) 



El = - { ——Hoo + Hii + H22 + H33 ] E2 



2 \ h J 2v 

^ Hn + H22 + Hs3 ^ -Hn-H22 + 3cos2e{~Hn + H22) + 2H^3-6sm2eHu 
Es = E, = 

(24) 
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us summarize the salient features of these equations: 

• The 15 equations £^'^ = spht up into five sets, decoupled from one another. 



The B and C sets fjl6|) -()20 |) involve only the Z2-odd components of the metric, and so 
it is inevitable that they are homogeneous. We may set Bi = B2 = C = and focus 
on the A, D, and E equations. 

The A equation (fT3j) happens to be identical to the dilaton equation of motion up to 
a factor multiplying the source term, so we may borrow directly from ^21 to find its 
solution. The C equation is also the same as the dilaton equation except that it 
is homogeneous. 



The B and D equations, p7|) - (fTS|) and (P^ - (f^^ . are identical except that the former 
are homogeneous and the latter are not. Each set involves one constraint and two 
second order equations of motion. 

The E set (j24j) - (j26j) involves four second order equations of motion and three con- 
straints. 

The total momentum K enters the equations of motion only as a multiplicative factor 
on the non-derivative coefficient matrices and through the source terms. Elaborations 
of the WKB method may therefore be suitable for approximately solving the equations 
at large K, and series solutions in K may be used at small K. Section 13.51 includes 
further discussion of small K approximations. 
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3 Analytic approximations 



Although the discussion at the end of the previous section allows us to set the three odd 
Hmn to 0, and correspondingly i?i = i?2 = C = 0, we sometimes refrain from doing so in the 
following discussion of limiting forms of the equations. 

3.1 Near the boundary 

By solving the equations of motion in a series expansion in y one obtains the leading forms 
X = -^y^ + Qfy^ + Rx X = A, B^, C, D,, D^, E,, E^, ^3, ^4 , (27) 

where 

PA = Pry^ = Po^=PE, = -PE, = l Pe,=1+V^ Pe, =v\3cOsH-1) 

(28) 

Pbi = = -Pc = . 

The Q^* are integration constants related to the VEV's of the stress tensor. The Rx are 
integration constants which can be set to zero because non-zero values would correspond 
to deformations of the gauge theory lagrangian. The constraint equations imply relations 
among the Qx^: 

/^tot ^tot _ ^ ^tot _ ^tot _ n 



2iK 

:i-3v^cosH)Qf^+2Q'^l 



3iv{l + v^) cos 6* 
2K 



The meaning of the equations (j29p becomes clearer in terms of the original variables Hmn, 
whose series expansion near the boundary of AciS's-Schwarzschild includes the leading terms 

Hmn = -^y' + Q'^ly' + Rmn . (30) 

Pmn, Qmn) a^d Rmn are linear combinations, respectively, of Px, Q^x", and Rx, as can be 
deduced by inverting the relations (jHl), (fTHjl . (fT!I|) . (j2H), and In particular, after setting 
Z2-odd quantities to zero and using to eliminate Q^^l, QeI, Q^eI^ QeI favor of 
g7, Q'£l = Q^^\ and Q'°l = Q^°\ one obtains 

Qmn OimnQ A ~^ ^mnQ £) ~l~ ^ninQ E ~^ Pmn 

(31) 
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where 
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(33) 



/ -4 

4-y cos^ 6* 
2?; sin 29 

\ 



4?; cos^ 9 2v sin 26* 

4eii (1 - 3r;2cos2 6l)sin26l 
3t;2cos2 6l)sin26l 4e22 

-2 + 2v'^ cos"^ 9 



en = 1 [-1 + (1 + 1,2) cos^ - 3v^ cos^ 0] 
622 = ^ cos^ 6'(-l - 2't;2 + cos^ (9) 
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2v tan 9 



-3 + t;2 + (1 - 3z;2) ^os^ 6* [-2 + (1 - 3i;2) cos^ 
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The Qj^Jii are integration constants which are proportional to entries of (Z^^): 



\ mnl 



Ql 



tot 
mn ' 



(35) 



(36) 



What we mean by {T^^ is a co-moving Fourier coefficient of the quark's contribution to 
(Tmn)- The overall factor in (jHEI) can be determined through ffist principles along the lines 
of [21], but we find it more instructive to obtain it heuristically by considering what may 
seem at ffist to be a digression: Ad^s-Schwarzschild in axial gauge. 
Defining a new radial variable q through 



y 



gV4 



(37) 
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one finds that the hne element Q becomes 

ds^ = -^^(-dt^ + dx^ + zjrdq'^) + h^^dx^dx'' 

. r3-gV4 ^ ^^^^ 

V=4;^diag 1^^-^,1,1, 1,0 

This is indeed an axial gauge description of AtiS's-Schwarzschild because h^o = 0; note 
however that (j38p is an exact rewriting of 0- 

On general grounds, the stress tensor of the boundary theory must be proportional to 
the coefficient of q'^ in hmn- But in the case of AdS's-Schwarzschild, there is a pre-existing 
determination of the energy density and pressure based on |25j : 

l=P= y iV^T^ . (39) 

Therefore we conclude that 

(T^„,) = ^N'T' hm -^^,iq'h,^r^) • (40) 

Returning to the setup with a string dangling into A(i5'5-Schwarzschild means that on top 
of the "perturbation" /i^,^ that deforms AdS^ into ArfS's-Schwarzschild we must add an 
additional perturbation, namely the h^j^^y whose Fourier coefficients are given in ()11|). The 
result po|) applies unchanged, except that the limit exists only after certain divergent delta- 
function contributions have been subtracted. After using (p?n|) and the standard relations 
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N\l = 4n'L' ^ = ^gl,,N (41) 
one obtains 

(r_) = ^iV^T^ diag{3, 1, 1, 1} + 1 ^(T^Je4^^(^^-*)+^-^+^3x3]/., ^42) 

where (T^„) is indeed given by 

Because Tmn is conserved and traceless (the latter due to conformal invariance), one 
expects that K'^Q^^^ = where 

K"" = (vKi Ki o) , (43) 
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and also 

tr Q*°* = -g*°* + Q\f + Ql°^ + g*3°3* = . (44) 
The tracelessness condition is indeed satisfied, but conservation fails: instead, 

K^'QTn) = f -1 o) . (45) 

The result fj45p is independent of Qa, Qd, and Qe- that is, only the last term in the 
decomposition (jHT|) fails to be conserved. 

The non-conservation (|45p could have been anticipated. It is the manifestation of the 
energy-momentum imparted by the quark to the thermal medium. The quark is prescribed 
to travel with constant velocity, so it does not slow down as it loses energy-momentum. 
The non-conservation (^3)) should precisely reflect the external force required to keep the 
quark's momentum from changing, which is minus the drag force (PQ). This argument is 
formal because the quark's mass is infinite, hence so is its momentum. But changes in the 
momentum, and therefore forces, can be finite. To verify that (0) can be recovered from 
fl45|) . consider some finite region V of R^. The external force on this region is 

F^ = ^ [ d^x{T^^)+(f d^ani{T^)= [ d^xd„,{T'^^) . (46) 
dt Jv Jdv Jv 

Here and in the following, i and j are three-dimensional spatial indices, while m and n are 
3 -|- 1-dimensional Lorentz indices. The first term in the middle expression of (j46p is the 
rate of change of momentum in this region, and the second term is the rate of escape of 
momentum through its boundaries. Using ()42|) . one obtains 



= [ d^X [ ^ J^^(rprnj^i[Ki{x^-vt)+K2X^+K3X^]/zH 

Jv J (271)^ zh ^ (47) 
i^m = [-vK^ K2 Kg) • 

The expression for in ()47|1 is equivalent to except that we have not specialized to 
K2 = K± > and = 0. Now take the limit where V covers all of so as to obtain the 
force on the whole system. Performing the x integral first, one obtains 

= izl I d'KK^{T^^)e-^^''^'6{K,)6{K,)6{Ks) = izl lim K„,{T^^) (48) 

J i?-»0 

where we have anticipated that K.m{T'^'') has a smooth limit a.s K Indeed, using (jSHI) 
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and then ()45|) we arrive at 

Vl — i?-»o 2 ^i_^2 

As promised, this is minus the drag force (0). 

To restore conservation of Tmn, one may add to it a "counterterm:" 



where, after passing to a co-moving Fourier description. 



K-r^n = -k-{t:^J = -10 0). (51) 



A solution to which is also traceless is 



r^lf Tif i 2 (l-i;2)3/2 1 _2t; 1 + 1;^ 



with other components of Tmn vanishing. Using 
d^K 1 



(27r)3 J'STi 



one finds 



(52) 



It would be cleaner if Tmn had delta function support at the location of the quark, but this 
does not appear to be possible: 7^ would then be analytic in Ki and K±, and there are 
no analytic solutions to (j5ip . The form (j54j) of the counterterm indicates an unphysical 
"string," wholly in the boundary theory, that pulls forward on the quark to counteract the 
drag force. 

The upshot of this somewhat extended discussion is that the original non-conserved form 
(|3ip captures the dynamics of dissipation and is non-conserved because it leaves out the 
external motive force that keeps the momentum of the quark from decreasing. 
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3.2 Near the horizon 

Near the horizon of Ac/S's-Schwarzschild, and for vKi ^ 0, the leading approximations to 
solutions to the equations of motion ()15|) and (j^^ are 



+ (^^^ yi-y) + Un (^-^^'-y^'j (1 - ,)--^/^ (56) 



16 + (vKi)^ l-ij^X ^1 1 + i^K- 



(57) 
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\ -^K\l-y) J 



+ Ve{1 - yT"'" 





VV 



(58) 

In these solutions, Ux-, Vx, and T^^ are integration constants. Near- horizon solutions to the 
B and C equations are identical, respectively, to the D and A solutions (j^Uj) and (jSH)), except 
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that the particular solutions are zero. 

For each set of equations, the solution multiplied by Ux is infalling (meaning that gravi- 
tons are falling into the black hole), while the solution multiplied by Vx is outfalling. The 
solutions multiplied by T^*^ are neither infalling nor outfalling but can be categorized by 
their regularity properties at the horizon. The standard boundary condition imposed at a 
black hole horizon is that outfalling modes must vanish: Vx = for X = A, B, C, D, and 
E. 

The constraint equations and imply 

^^2) ^ y|) ^ y(4) ^ y(5) ^ Q _ (59) 

The solutions in ()56|) and (j3Hj) multiplied by Tj^\ Tj^\ and T^^ are in fact exact solutions 
to the equations of motion for all y. Note that the exact solutions do not overlap with the 
ones removed by (j59|) . This suggests that in the coupled systems of equations of motion and 
constraints for Di and Ei, it may be possible to make further reductions of order. We have 
not pursued this avenue, but it might facilitate future numerical studies. 

To understand the boundary value problem that determines (T^„), it is useful first to 
review the counting of integration constants, constraints, and boundary conditions: 

• The ten second order equations of motion have 20 constants of integration which must 
be fixed in order to specify a unique solution. 

• Ten constants of integration are fixed by requiring Rx = at the boundary (no defor- 
mations of the gauge theory). The other ten are the Q^x^, which are linear combinations 
of entries of (T^„). 

• Five relations among the follow from imposing the constraints at the boundary. 

• Five more boundary conditions must be imposed at the horizon to suppress the out- 
falling solutions. 

Evidently, the number of constraints plus boundary conditions equals the number of inte- 
gration constants in the equations of motion. So the boundary value problem is well posed. 
All the integration constants are complex, and the constraints and boundary conditions are 
too. 

Similar counting of integration constants can be done after dividing the equations into 
the decoupled sets, A through E. Let us include this counting in a summary of the numerical 
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algorithm we used. In the A set, one must impose Ra = 0. If we supply in addition an ad 
hoc value q for Q^*, then Cauchy data has been specified at the boundary. More precisely, 
approximate Cauchy data can be specified at a finite but small value y = Ho hj setting 
Ra = and Q^^^ = q equal to its ad hoc value and using (jTfj) to determine A^i/q) and A'^uq). 
The second order equation of motion ()15|) can then be integrated numerically to y = yi, 
where yi is close to 1. The numerical solution can then be fit to the asymptotic form (j55p . 
and values of Ua and Va can be extracted. Because all the equations are affine (meaning 
linear with inhomogeneous terms), Va (as well as Ua) is an affine function VA{q) of the ad hoc 
value we supplied for Q^'- The equation VA{q) = may easily be solved for the physically 
meaningful value of Q^*. 

For the D and E sets, the situation is only slightly more complicated. After setting 
Rx = (see ()27|) ) and imposing the constraint (for D) or constraints (for E), there is only 
one degree of freedom left at the boundary, which we can fix by supplying an ad hoc value for 
the quantity or Q^* that enters ()H Ij) . Cauchy data for the equations of motion can be 
generated aX y = yo, and after numerically solving the equations of motion, the integration 
constants t|\ Vx, and Q^* can be extracted by matching numerics to horizon asymptotics 
at y = yi. To determine Q^* or Q^*, one solves an affine equation Vx{q) = 0. It would have 
been numerically more efficient to eschew one of the equations of motion in the D set and 
all but one in the E set in favor of the first order constraint equations. But we found it a 
useful check of numerical accuracy to evaluate at y = yi the T^^ which are required by (jSHI) 
to vanish. 

The method of obtaining an affine function at the horizon by first specifying Cauchy data 
at the boundary was described in (211 for graviton perturbations in ArfS's-Schwarzschild in 
the absence of the trailing string. 

3.3 Large K behavior 

The large K behavior of Q^^n is dominated by the near-field of the quark. For v = 0, 
this field is entirely color-electric, and it is perfectly Coulombic because of the conformal 
symmetry of A/" = 4 super- Yang-Mills theory. These observations alone, together with the 
radial symmetry and the conservation and tracelessness conditions, fix the v = form of 
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Q^mn overall prefactor: 



y^near 



IT 



24\K\ 



/2K^ 
K| + Ki -K1K2 
-K1K2 Kl + K\ 
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(60) 



In section 123 we explain how to fix the overall prefactor in (jHOJ. 

Q'^n for f 7^ can be obtained by applying a Lorentz boost to (plj) . After this is done, 
one may define 



Qmn Q 



tot 
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near 

mn ' 



Then — > as i^T — > 00.^ More useful in section HI will be the equivalent forms 



(61) 



Qx = QT - QT" x = A,D,E 



(62) 
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24 VI - cos2 9 



(63) 



Recall that Qd = and Qe = Qei- To derive (jU^ one must compare the conserved 
terms in (|!?T|l with the Lorentz-boosted version of ()6()|1 . with = and Jir2 = K± in the 
rest frame of the thermal plasma. 



3.4 Normalizing the near field 

The prefactor can be fixed by observing that the equations pop for v = 0, h = 1, and 
^ = 0, as appropriate for a static string in AdS^, can be solved for K_i_ = by setting all 
the Hmn = except for Hqq, Hu = 2y^/9, and H22 = -f^ss = = -2//3, where f{y) 

satisfies 

f = y- (64) 

This is precisely the equation satisfied by the K± = Fourier modes (pK of the dilaton 
sourced by the same static string configuration: see (17) of ^21- Although (jMj) is non-trivial 

^Actually, can be arranged to have an arbitrary A'-independcnt limit for large K by adjusting the 
choice of particular solution. This corresponds to adjusting a subtraction scheme for the infinite self-energy 
of the external quark. The form of solutions specified in (|30|1 . where the particular solution is assumed not 
to have a quartic term, is a sort of holographic minimal subtraction scheme. 
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to solve directly, position space methods are available to extract the dilaton profile |27j . 
From them one can Fourier transform back to find = 7r|A'i|/16, where is defined 
through the asymptotic behavior 

+ (65) 

From Hqo = — 4//3 it follows that Q^q'^'^ = — 4i?i^/3; hence the prefactor in (j60|) . Comparing 
(23) of 1121 (EEI) of the current paper, one may conclude that 

{Too) = 1{Of2). (66) 

This is a positive quantity because Op^ ^ ti E"^ for the static quark. We do not know how 
to account for the factor of 4/3 in (jM]) . 

It is instructive to examine the same static quark solution using the equations (fT^ - ()26|l . 
Taking h = 1 and ^ = in these equations poses no difficulties. K± = means 6 = 0, which 
appears to lead to difficulties in ()14|) and (PT|) (the definitions of A, Di, and D2 in terms 
of Hmn)- But the inverse relations expressing Hmn in terms of the ABODE variables are 
entirely non- singular in the limit i^x 0: they read 

2 2 

Hqq = "g-^i + g-^a Hqi = 2vE2 

2 2 

H02 = Hii = -Es — -E4 

2 1 ^ ^ 

H12 = H22 = —E3 + -E4 



3 



= g-E's + -^E^ 



We have omitted expressions for the Z2-odd components of the metric in terms of Bi and 
C because all these quantities vanish in the solution we're interested in. The equations of 
motion and constraints for A, Di, and Ei are also non-singular in the limit K± 0. The 
equations of motion and constraints for A and Ei are non-singular if one additionally takes 
V 0, but the constraint for Di is not. Therefore in (j68p - (j72j) we partially quote and 
partially solve the equations of motion and constraints for A and Ei after having set v = 0, 
but for Di we keep v finite. 

dl-hy-K', A = y (68) 
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takes precisely the same form as ()64|) . 



W 



AKi 1 - ^2 



D2 = d, 



W 



is the general solution of the D constraint consistent with the requirement that Di 
y ^ 0. The function di satisfies 



di = y; 



which takes the same form as ()64|) except for the replacement Ki — > Ki^Jl — v"^. 

2 



w 4 

-c/i = -ei H 

3 9 



2 y-^ 

En = —E-] = -ei Ea = — ei 

^39 3 



9 



is the general solution of the E constraint consistent with the requirement that Ei 
y ^ 0. The function ei satisfies 



9y - -dy - Ki 



ei = y: 



(69) 



as 



(70) 



(71) 
as 

(72) 



which again is precisely the same form as (jMjl . 

Because (j67|) involves only the Ei, it was superfluous to explicitly solve the A and Di 
equations in (jUHj) . But it is a worthwhile check to ensure that the quantities Di, though 
singular in the limit v ^ (as well as in the limit Ki — > 0), cause no problems for Hmn- 
Indeed one recovers the results for Hmn stated briefly around (j64|) . 



3.5 Small K behavior 

A series expansion in small K allows some progress to be made on solving the equations of 
motion. We will first focus on the simplest case, namely the A equation. Plugging 



A = ao + Kai + K^a2 + ■ ■ ■ 



(73) 
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into p5|). one can find the differential equations satisfied by each aj. The first three are: 

'4dy\d,a, = (74) 
h y-^ h zh 

h y f cos cos^ Q — h 
= _ 

Evidently, these equations are solvable through repeated integration. Integrating is easy, 
and after matching to the boundary asymptotics with Ra = and the horizon asymptotics 
with Va = 0, (both suitably expanded in K) one obtains 

ao = -l log(l + y)-] log(l +y')+'-r log ^ = + ^ + 0^) , (75) 
2 4 4 1 + ty 6 4 

which implies = 1/4 + 0{K). Higher order corrections to can be found by solving 
the corresponding differential equations for each of the aj's, and matching the solutions to 
the horizon and boundary asymptotics. To order OlK"^) we have obtained 



1 ?■ loe- 2 

QT = ^ - -^vK, + ^{Gn - 121og2) 



2 , 2^2-12(l0g2) , 2 



(76) 



l-v' + tf ^ , ' cos^ 6 

67r-121og2 



0{K'' 



A similar analysis can be carried out for the D and E sets. For the D set, for example, 
one first writes down series expansions in K for the second order differential equations 
and solves for the corrections to Di and D2 iteratively. By imposing the constraint 
and by matching, at each order, the small K solution to the boundary asymptotics with 
-Rdi = = s-iid to the horizon asymptotics with Vd = 0, one can then solve for the four 
integration constants and the corresponding corrections to Q^*, Q^*, Ud, and Tj^^ . To 
linear order in K, we find 

i sec 6 sec^ 6 - 4f ^ iK sec^ 6 . ,^ , ^,5 2 ^ 

= - «r + ^4^^ [1 + P 2'" ™^ « (77) 

- (81og2)t;^cos^^l +0{K^). 
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Figure 1: The A(i5'5-Schwarzschild background is part of the near-extremal D3-brane, which 
encodes a thermal state of A/" = 4 supersymmetric gauge theory j2S]- The external quark 
trails a string into the five- dimensional bulk, representing color fields sourced by the quark's 
fundamental charge and interacting with the thermal medium. 



Similarly, for the E set, we find 
3iv(l + v'^)cos9 



3t;2 cos2^[2 + t;2(i _3cos2^)] 



2K {I - ?>v'^ cos^ 6) 
3u^(l + w^)cos6' 



2fl 



"iv"^ COS^ I 

1 



0{K) 



(7J 



2K 



;i -3t;2cos2^) (1 - ^^S^ ) -ivKcosO 



0{K). 



The striking feature of the expression is the singular behavior at ^ = cos~^(l/t'v^), which 
is the Mach angle. From this we may conclude that there is a sonic boom in the thermal 
medium involving large amplitude but small momentum fields. In the second expression, we 
have shown how the 0(1) term may be "resummed" into the leading 0{1/K) expression so 
as to blunt the singularity into a form resembling a Lorentzian lineshape. 



4 Results of numerics 

Let us briefly recap the five-dimensional gravitational calculation that has been our main 
focus so far. The trailing string of |2] sources the graviton, which propagates classically 
in AdS's-Schwarzschild with purely infalling boundary conditions at the black hole horizon. 
The graviton's behavior near the boundary of ArfS's-Schwarzschild determines (Tm„) in the 
boundary gauge theory. Thus (Tmn) is a shadow (other authors might prefer the term 
"hologram" ) of the trailing string. See figure ^ 
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a) Re Qa(Ki,Kp) for v=0.95 b) -Im Qa(Ki,Kp) for v = . 95 c) Kp |Qa(Ki,Kp)| for v=0 . 95 




Figure 2: Contour plots of Q^, Q^, and for v = 0.95. The darker regions are more 
positive. All components of (T^„) can be deduced from Q^, Q^, and using (|HT|) . 
and dnS). All three go to zero at large K. The momentum vector K = k/nT can be read 
in GeV/c if one chooses T = 318 MeV: see (j79|) . The range of momenta in each plot was 
chosen to show the most distinctive structures. The region outlined in gold in (i) is plotted 
in more detail in figure Et- 
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a) Kp I Qe (Ki, Kp) I forv = 0.75 b) Kp | Qe (Ki , Kp ) | forv = 0.9 




01234 01234 



Kl Ki 

Figure 3: Contour plots of -ft'j.lQ^I for various values of v. is proportional to the i^'-th 
Fourier component of the energy density after a near- field subtraction: see (jHTj) . (jHUj) . and 
(jU^ . The phase space factor K± arises in Fourier transforming back to position space. The 
green line shows the Mach angle. The red curve shows where K±\Q^\ is maximized for fixed 
K = The blue curves show where -R'.lIQ^I takes on half its maximum value for 

fixed K. 

Our aim is to describe {Tmn) in the boundary theory. We will focus on Fourier coefficients 
for X = A, and E. As reviewed at the end of section Q these quantities are Fourier 
coefficients of linear combinations of entries of (Tmn) with a near-field subtraction. Our 
numerical algorithm is outlined at the end of section 1X21 It was implemented primarily using 
Mathematica's NDSolve. To achieve good accuracy, it was necessary to develop asymptotic 
power series solutions to a considerably higher order than shown in sections 13.11 and 13.21 
Experience as well as common sense suggest that large K regions become more numerically 
challenging. We believe we have adequately met this challenge, partly by allowing arbitrarily 
many steps in NDSolve and calculating with a working precision of 30 digits (i.e. roughly 
twice the standard double precision of modern PC's). Another challenging region is small 
Kl, where the outfalling and infalling solutions are nearly constant until y is very close to 
1. Experience suggests that at most a narrow region with Ki -C is problematic. 

We found excellent agreement between the numerically computed and the analytical 
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a) Kp|QH:(K:,Kp) I f or v=0 . 4 b) Kp | Qe ( Ki , Kp ) | forv=0.5 




0.02 0.04 0.06 0.08 0.02 0.04 0.06 0.08 



Ki Ki 

Figure 4: Contour plots of K^\Q^\ for various values of v at low momenta. The green 
line shows the Mach angle. The red curve shows where /T^IQ^^I is maximized for fixed 
K = Kf + Kj_. The blue curves show where K±\Q^\ takes on half its maximum value for 
fixed K. 
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approximations Q^^^ for large K and for X = A, D, and E. Only the subtracted quantities 
Qx = Q^x ~ Qj{^^ appear in the plots in figures 121121 For and Q^, we also found 
excellent agreement with the small K asymptotics (fTUI) and (fTTj) in the expected ranges. For 
Q^, the nearly singular behavior near the origin is difficult for numerics to capture. This 
difficulty shows up in the ragged contours in figure |3t,d. The problem is not numerical 
error in evaluations of at individual points; rather, the ragged contours in figure lltjd are 
due mostly to imperfect interpolations over a grid of limited resolution. Indeed, individual 
evaluations of Qf for K = ^ Kl + K\ = 0.08 agree with the small K asymptotics (f7S|) at the 
level of about a percent. A high-resolution plot of -R'.lIQ^I at = 0.08 is shown in figureEfe'- 
In this plot, the results of numerics are visually indistinguishable from the analytic form ()78|1 . 
Even in the coarser-grained numerical evaluations of shown in figure |31 agreement with 
(f7H|) was good a distance 6K ~ 0.015 away from the central ridge. 

Two qualitative features visible in figures 121121 are worthy of note. The first is the high 
momentum ridges, which are most distinctive in figure |2t. This is the same feature that was 
noted in jT21; indeed, of this paper is identical to Bk of ^21- High momentum ridges are 
also present in and Q^. For v = 0.95 and v = 0.99, we find empirically that ~ 
on the high momentum ridges. A more approximate relation for a similar region of momenta 

is Qf = sgf . 

The second feature worthy of note is the sharp structures at low momentum in fig- 
ures |2|;;h,i. A more detailed view of these structures is shown in figure El As we will discuss 
in section the lobes in figure El are suggestive of high angle emission of particles in en- 
ergy ranges accessible to experiments at RHIC. The lobes become narrower as one passes 
to small K, corresponding to momenta much less than the temperature: see figure |31:;,d. 
Low momentum is the hydrodynamic limit, so it is gratifying to see a highly directional 
feature corresponding to a sonic boom. Figures lll:;,d thus serve as visual confirmation of 
the appearance of a sonic boom that we anticipated based on (f7S|). Figures |l^,b show what 
happens when the velocity of the quark falls below the speed of sound ~ 0.577 in 

the thermal medium. Evidently, there is still directional emission, but it becomes abruptly 
less focused. The peak amplitude also decreases abruptly. Intriguingly, the drag force ((H) 
behaves completely smoothly as one passes through v = 1/ ^/S. 

Readers wishing to examine our results more quantitatively are referred to [2H1- 
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5 Discussion 



Because all the calculations in this paper were carried out in the framework of five- dimensional 
supergravity coupled to a classical string, the gauge theory results are accurate only to lead- 
ing order in large and large Qym^- Large 't Hooft coupling is largely inaccessible to 
standard techniques of finite-temperature quantum field theory, with the important excep- 
tion of lattice gauge theory. But finite-temperature lattice methods are not well-adapted to 
real-time dissipative phenomena, in contrast to AdS/CFT, which provides ready access to 
both static and dissipative properties. Moreover, the AdS / CFT prescription for computing 
gauge theory observables is conceptually the same at all energy scales, giving some advan- 
tage over hydrodynamical approximations that are best justified in the infrared limit. So 
AdS/CFT occupies a unique niche in the range of tools available for understanding strongly 
coupled gauge theories at finite temperature. Its principal drawback, at least as we use it 
in this paper, is that the dual gauge theory is SU{N) A/" = 4 super- Yang-Mills, which in 
some ways is quite distant from real-world QCD. Within the limitations that we have de- 
scribed, the calculation of (Tmn) provides a fairly comprehensive description of dissipation 
from the heavy quark. All possible gauge interactions are included, in particular secondary 
interactions with the thermal medium of energetic particles radiated from the quark. 
For the sake of definiteness, let us set 

T = -GeV = 318MeV. (79) 

TT 

We understand this number to be in the upper range of temperatures for the quark-gluon 
plasma (QGP) produced at RHIC. It is a convenient choice for us because the Ki and K± 
axes in figures |21 El and El can then be read in units of GeV/c. 

In ^2] we suggested that the high momentum ridges might be evidence that the strongly 
coupled thermal medium enhances fragmentation near the kinematic limit. But one of our 
warnings was that one should compute (Tmn) before making definitive statements. In light 
of the pronounced directional peak in i^xIQ^I inclined to regard the 

high momentum ridges as less immediately important to attempts to compare string theory 
calculations to recent experimental results on jet-quenching. It is plausible that these ridges 
are the expression in Fourier space of a sharp "prow" of color fields supported near the quark. 
Moreover, we must bear in mind that at the typical energy scale ^ 10 GeV where the high 
momentum ridges are pronounced, QCD is no longer particularly strongly coupled, so there 
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is less justification for a connection witli tlie supergravity approximation in AdS/CFT.^ 

Our most striking results are the directional lobes in Kx\Q^\, as seen in figure |21 Recall 
that the factor of K_<^ is appropriate because it is the measure factor arising in a Fourier 
transform back to position space after the azimuthal integral is performed. Recall also 
that is directly proportional to the K-th Fourier coefficient (Tqq) of the energy density 
with the Coulombic near field subtracted away, whereas and are combinations of 
components of (T^„) and non-conserved terms: see (jH!)), (jHBj) . and (jU^ . From figure E] we 
conclude that in strongly coupled A/" = 4 gauge theory at finite temperature, directional 
emission from a hard probe is present, but not sharply focused, at momenta several times 
the temperature. This seems to us an intuitively appealing conclusion: Rescattering effects 
broaden the directionality of the "wake" in Fourier space. 

In figure El we show a small sampling of our numerical results in a format more suggestive 
of a comparison with experimental results [221 122| on the splitting of the away side peak 
in di-jet hadron pair correlations. What we find attractive is that at momenta comparable 
to the window 1 GeV/c < Pt < 2.5GeV/c of transverse momenta for the partner hadrons, 
there are broad peaks in the Fourier components (TJ^) of the energy density: see figure int,d. 
These peaks are not unlike the ones seen at A0 ^ 2 in gold-on-gold collisions more central 
than 60%: see figure 2 of j2Sl- They are quite different from the narrow peak that we 
find at 80MeV/c (figure Et^). At 80MeV/c, which is about a quarter the temperature, a 
hydrodynamic description is probably justified, and we should interpret the narrow peak as 
a sonic boom. The analytic form 

,j.K. ^ 7r^T^V^h^3vil + v')cose 1 

^00/ 2iK (i_3^,2cos2^) (1- Hg^) -^t;i^cos^ ^ 

(80) 

which follows from (jSEI), and (ffH|) . is highly accurate in the infrared limit. One can see 

from figureElthat for v = 0.95, ()80|) loses validity around K ^ 1, corresponding to 1 GeV/c. 

In the interesting region of 1 to a few GeV/c, (T(^) decreases significantly more quickly with 

increasing K than the approximation ()80j) would indicate. This falloff may be a positive 

feature in comparing to data. 

Let us enumerate the reasons to treat with caution the connection we allege between our 

•^When hadron pair correlators are plotted with higher momenta windows for the hadrons, we understand 
from an experimental colleague that the away side peak reappears |29) . The high momentum ridges might 
be relevant to such correlators: forward emission is indeed what they imply. But it is perhaps more plausible 
to attribute the reappearance of the away-side peak to away-side partons that have enough energy to punch 
through the QGP with only modest deflection. 
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a) Kp|QE(Ki,Kp) I for K = 0.08 b) Kp|QE(Ki,Kp) | for K = 0.8 




Figure 5: -R'xIQ^I at fixed K = ^ Kl + K\ as a function of angle, for v = 0.95 and for 
various values of K. To facilitate comparison with di-jet hadron pair correlations, we have 
parameterized the angle as A0 = ir — O, where 6 = tan~^ Kj_/Ki. With the usual assignment 
T = 318 MeV (see ^7^ ). K can be read in units of GeV/c. In each plot, the solid curve is 
from numerics; the dashed curve is the analytical approximation (j78|) : the green line indicates 
the Mach angle; the red dot is at the maximum of Kj^\Q^\\ and the blue dots indicate the 
points where K^IQeI is half of its peak value. 
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AdS/CFT calculations and the RHIC results on away side jet splitting. 

1. is the hardest quantity to compute of the three that we investigated: the system of 
equations is more formidable, and the nearly singular behavior of near the origin 
makes numerical evaluations less stable. Moreover, our analytic approximation (j78|) 
is not as precise as for and Q^, meaning that we have less extensive checks on 
numerics. 

2. The A(j) in figure 2 of [23| is the separation in azimuthal angle, whereas in our figure El 
it is TT minus the angle between the emission direction and the motion of the heavy 
quark. 

3. The broad peaks at A0 ^ 2 in are distinctive only after a subtraction related to 
elliptic flow. 

4. The peaks of di-jet hadron pair correlations are closer to A0 ^ 2 than to the peak 
angle A0 ^ 2.4 in figure El 

5. The experimental studies fl^ (221 do not include heavy quark tagging, so most of 
the away side partons are presumably light quarks or gluons. But perhaps, for high- 
angle emission, what matters most is not the quark mass but simply the color current 
associated with a hard parton. 

6. After a parton leaves the QGP, it fragments, and then its fragmentation products 
must be detected. We do not have the expertise to add these important aspects of the 
physics to our calculations. 

7. The QGP cools, expands, and hadronizes, and its equation of state changes with time as 
a result. The conformal result Cs = 1/ -\/3 is likely to be a reasonable approximation only 
in the QGP regime, at temperatures significantly above the deconfinement transition. 
As remarked in jHOj, a steeper emission angle results from a time-averaged speed of 
sound that could be as low as Cg ~ 0.33. It might be possible to partially mimic the 
changing equation of state by some deformation of AdSs-Schwarzschild, but it's not 
clear that the result would have the same status as a first-principles calculation that 
can be claimed for our analysis. 

8. One of our many idealizations of the true experimental setup is that we replaced the 
QGP by a thermal medium of infinite extent. This could mean that we are exaggerating 
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the effects of secondary rescatterings. It would be desirable to have some position space 
representations of components of the stress tensor to address this issue. 

9. It remains to us a deep mystery when and why strongly coupled A/" = 4 gauge theory 
should be directly compared with real-world QCD. Doing so somewhat above the 
deconfinement and chiral symmetry breaking transitions is clearly the best hope. But 
we return to the basic conundrum: are near-extremal D3-branes merely an analogous 
system to the QGP, or can they capture the dynamics sufficient precisely to be a useful 
guide to RHIC physics? 

It is clear from figurelHlthat as one passes to higher momenta, the peak emission direction 
becomes more forward, although at the same time the peak keeps broadening. It would be 
interesting to compare the dependence of {T^) on both the magnitude and angle of K with 
two-dimensional histograms of Pt and A0 for partner hadrons. 

As K increases beyond the range shown in figureEl one eventually passes into the region of 
high momentum ridges, which become more and more forward as v 1.^ As we understand 
the experimental situation, the away side peak reappears as one increases the momentum 
window for the hadrons. An optimistic read of this situation is that AdS/CFT calculations 
may have some relevance up to an unexpectedly high range of momenta; but perhaps it 
is more reasonable simply to suppose that sufficiently high-momentum partons can punch 
through the QGP without much deflection. 

There have been other notable theoretical efforts to understand the splitting of the away 

side jet. An account of the sonic boom picture can be found in [201 • Investigations of the 

Cerenkov radiation have been pursued in I33j. The conical flow picture has been seen in 

[SI] using the hydrodynamical evolution of QGP, and in [33] using linear response theory. In 

[53] the Mach cone picture appears only in the strongly coupled QGP. In comparing with 

these more phenomenological works, it must be admitted that we have made dramatic and 

risky idealizations of the experimental setup. Yet, despite the potential stumbling blocks, 

it is exciting to see a simple type IIB string theory construction approaching quantitative 

comparisons with a data-rich experimental field. 

^We thank J. Casalderrey-Solana for pointing out to us that the peak angle of the high momentum ridges 
decreases roughly as I/7, and for the interesting remark that this behavior may signal some connection with 
the Landau-Pomeranchuck-Migdal effect (see for example 
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